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ACCEPTANCE  CONTROL  CHARTS 
BASED  ON  THE  EXACT  AND  NORMAL  APPROXIMATIONS  TO 
THE  BINOMIAL  DISTRIBUTION 


by 

Carlos  Amado 
Richard  S.  Leavenworth 
Richard  L,  Scheaffer 


ABSTRACT 

Procedures  are  developed  for  finding  the  sample  size  and  control  limit 
for  Acceptance  Control  Charts  for  proportion  of  nonconforming  units  using  the 
exact  binomial  distribution,  the  standard  normal  approximation  to  the 
binomial,  and  a  normalized  arcsin  transformation  of  the  data.  The  user  must 
select  an  Acceptable  Process  Level  and  a  Rejectable  Process  Level  and  the 
associated  risks  for  each.  The  approximation  methods  are  compared  to  the 
exact  method  over  a  wide  range  of  design  specifications.  It  was  found  that 
the  arcsin  transformation  Is  considerably  more  accurate  than  the  standard 
normal  approximation  and,  although  more  complex  to  figure.  Is  preferable  If 
the  user  Is  familiar  with  small  scientific  pocket  calculators.  If  a  micro¬ 
processor  or  minicomputer  Is  available,  the  exact  binomial  may  be  used  with 
ease  to  achieve  at  least  the  stipulated  risk  protection  desired.  FORTRAN 
programs  for  the  three  formulations  are  Included. 


ACCEPTANCE  CONTROL  CHARTS 
BASED  ON  THE  EXACT  AND  NORMAL  APPROXIMATIONS  TO 
THE  BINOMIAL  DISTRIBUTION 


INTRODUCTION 

The  control  chart  for  fraction  rejected,  or  p-chart,  probably  is  the  most 
widely  used  of  all  control  chart  procedures.  It  may  be  applied  to  one  or  more 
than  one  quality  characteristic  whether  measured  on  a  go,  not-go,  basis  or  as 
variables  measurements.  So  long  as  the  result  of  an  inspection  is  to  classify 
an  item  as  meeting  specifications  (acceptable)  or  not  meeting  specifications 
(unacceptable),  a  single  p-chart  may  be  used. 

When  the  sample  subgroup  size  is  constant,  the  chart  for  may  be  used 
conveniently  since  it  records  the  actual  count  of  rejected  units  in  a  subgroup 
of  size  ^rather  than  the  proportion  rejected.  In  either  case,  the  binomial 
probability  density  function  may  be  used  to  model  the  process. 

The  standard  Shewhart  control  chart  for  places  the  upper  and  lower 
control  limits  at; 

n^  3  /  n'p(l-^ 


where : 

n£  is  the  average  count  of  rejected  units  for  a  series  of  subgroups 
of  constant  size,  n. 


/ np(l-p)  is  the  standard  deviation  of  the  binomial 


count  (a  ). 
'  n£ 


Thus  these  limits  are  standard  3-sigma  limits  of  the  Shewhart  control  chart. 
While  variation  beyond  these  limits  at  random  should  occur  very  rarely  indeed, 
the  same  probabilities  often  associated  with  random  variation  beyond  3  sigma 
limits  on  an  J  chart  should  not  be  applied.  Nevertheless,  exact  probabilities 


calculated  from  the  binomial  distribution  may  be  found  for  ^  control  charts 


with  3>sigma  limits  or  any  other  set  of  limits  which  are  a  multiple  of  o^. 


In  the  case  of  the  Acceptance  Control  Chart,  at  each  point  a  sample  is 
selected  a  decision  is  to  be  made  either  to  accept  the  hypothesis  that  the 
process  is  operating  as  specified,  at  the  acceptable  quality  level  or  better, 
or  that  the  actual  quality  level  is  beyond  an  acceptable  level  (a  higher 
value  of  £).  The  decision  criterion  is  the  upper  control  limit  on  np, 
i!B.  It  where  y  is  the  probability  of  acceptance  of  the  hypothesis 
with  respect  to  2.  using  a  normal  approximation  to  the  binomial  distribution. 


(Z  equals  3  for  the  standard  Shewhart  control  chart.) 
”Y 


This  paper  explores  the  development  of  Acceptance  Control  Charts  for 
binomial  counts  of  the  number  of  units  rejected.  Two  parameters  are  to  be 
found,  one  the  value  of  the  control  limit,  and  the  second,  the  appropriate 
(constant)  sample  size,  ju 

The  problem  is  formulated  in  three  ways,  using  the  exact  binomial  dis¬ 
tribution,  using  a  standard  normal  approximation  to  the  binomial,  and  using 
a  normalized  arcsin  transformation.  Analytical  results  of  application  of 
the  three  methods  are  compared  and  the  results  of  a  simulation  study  using 
computer-generated  synthetic  data  are  presented. 


PROBLEM  FORMULATION;  EXACT  BINOMIAL 


If  samples  are  being  drawn  from  a  continuous  process  generating  noncon¬ 
forming  items  at  a  constant  rate,  the  binomial  distribution  describes  this 
process  accurately.  The  probability  density  function  is 

^(f*!n,p)  =  (J)  p'^(l-p)"''^  (1) 


where 


r  =  no.  of  units  rejected 
n  =  sample  (or  subgroup)  size. 


The  probability  that  the  number  of  units  rejected  Is  less  than  or  equal  to 
some  fixed  value,  Is: 

P[r  <  c|n.p]  .  I  (")  p''(l-p)"-''  =  T  (2) 

r=0  ^ 

Limits  for  acceptance  control  charts,  like  acceptance  numbers  for 
sampling  plans,  usually  are  designed  by  selecting  two  points  on  an  Operating 
Characteristic  (OC)  curve,  as  Illustrated  In  Figure  1.  One  point  assures  that 
an  Acceptable  Process  Level  (APL),  pj,  has  a  high  probability  of  being 
accepted  of  at  least  1-a;  the  other  assures  that  a  Rejectable  Process  Level 
(RPL),  P2,  will  have  a  suitable  low  probability  of  acceptance  of  at  most 
a  and  p  are  the  design  risk  levels  associated  with  the  process  quality  levels 
Pi  and  P2,  respectively,  where  P2  Is  greater  than  pj. 

The  two  equations  needed  should  express  the  fact  that  It  Is  desired  that 
the  OC  curve  pass  through,  or  pass  as  close  as  possible  to,  the  two  points 
specified  as  the  design  criteria,  namely,  (£i»l-o)  and  These 

equations  are: 


P[r  <  c|n,p,]  =  1  (")  pJCl-Pi)""'*  >  1-a 
r=0 

(3) 

PCr  <  cln,P2]  =  (J)  PgCl-Pe)"”''  <  P 

(4) 

They  are  expressed  In  Inequality  form  because,  for  any  Integer  values  of 
^  and  £,  It  Is  unlikely  that  the  cumulative  probability  can  exactly  satisfy 
1-a  and  p.  As  the  result  there  are  an  Infinite  number  of  (n_,c_)  pairs  satis¬ 
fying  (3)  and  (4)  above  some  minimum  combination.  The  complete  formulation 
thus  requires  the  further  stipulations: 

minimize  n  and  c  1 


where: 


P2  >  Pi  and  1-a  >  p 


(5) 


Figure  1.  Design  Operating  Characteristic  (OC)  Curve  for  np  Chart. 


Referring  to  Fig.  1,  equation  (3)  assures  that  the  resulting  OC  curve  will 
pass  above  and  to  the  right  of  the  design  curve  at  the  point  1-a). 
Equation  (4)  assures  that  the  resulting  OC  curve  will  pass  below  and  to  the 
left  of  the  point  (£2»P)*  Stipulating  the  choice  of  the  minimum  (jt^.c)  pair 
assures  that  the  resulting  OC  curve  passes  as  closely  as  possible  to  the 
design  curve. 

The  value  of  of  course,  yields  the  constant  subgroup  size  to  be  used 

in  sampling.  Thus  the  ng_ chart  becomes  a  reasonable  and  easily  understandable 

« 

alternative  to  the  ^  chart.  The  value  of  ^1s  the  maximum  count  of  rejected 
units  that  should  lead  to  no  corrective  action  on  the  process.  Only  if  (c+1) 
or  more  units  are  rejected  should  action  be  taken.  However,  if  the  Acceptance 
Control  Limit  (ACL)  is  plotted  exactly  at  the  value  of  £,  the  user  may  become 
confused  as  to  whether  or  not  to  take  action.  In  accordance  with  the  rules  of 
control  chart  interpretation,  this  is  not  a  coin-flip  situation.  A  reasonable 
procedure  is  to  plot  the  Acceptance  Control  Limit  at: 

ACL  =  c  +  0.5  (6) 

thus  avoiding  any  confusion  in  chart  interpretation. 

PROBLEM  FORMULATION:  NORMAL  APPROXIMATIONS 

Working  with  the  binomial  formula  presents  a  number  of  mechanical  prob¬ 
lems  some  of  which  are  discussed  subsequently.  Suffice  it  to  say  at  this 
point  in  the  discussion  that  no  closed  form  solution  for  the  values  of  ^  and  £ 
exists.  Its  application  involves  repetitive  use  of  some  form  of  search 
algorithm.  Therefore  no  one  should  be  surprised  that  considerable  attention 
and  ingenuity  have  been  applied  in  the  development  of  useful  approximations 
to  the  cumulative  binomial. 


Johnson  and  Kotz  (1969)  provide  a  rather  extensive  survey  of  binomial 
approximation  techniques  and  Raff  (1956)  has  compared  the  accuracy  of  several 
of  them.  The  two  presented  and  compared  in  this  study  are  the  standard  normal 
approximation  and  a  normalized  arcsin  transformation.  The  standard  normal 
approximation  is  the  most  familiar  and  used  as  alluded  to  in  the  Introduc¬ 
tion.  It  is  easy  to  apply  since  ^  and  c_may  be  obtained  directly  with  the 
use  of  a  slide  rule  or  pocket  calculator  and  a  table  of  the  standard  normal 
curve.  The  arcsin  transformation  is  somewhat  more  complicated  but  easily 
adaptable  for  use  on  a  programmable  pocket  calculator. 

Standard  Normal  Approximation 

The  mean  and  standard  deviation  of  the  binomial  distribution  are: 

E(r)  =  np 
Op  =  /  np(l-p). 

The  distribution  of  the  standardized  binomial  variable 
Z  =  (r-np)/  /np(l-p) 

tends  to  the  standard  normal  distribution  as  becomes  large.  (See  Johnson 
and  Kotz,  1969.)  That  is,  for  any  real  number,  X: 

Lim  P[Z  <  X]  =  exp(-u^/2)  du 

n+“>  /  2ix 

which  values  may  be  found  in  a  table  of  the  normal  curve  or  solved  for  on  many 
programmable  pocket  calculators.  Thus,  if  given  the  value  of  an  ACL  (say, 
derived  from  the  exact  binomial,  non-integer,  and  equal  to  c^ +  0.5),  the 
probability  of  c_ or  less  occurrences  with  ^  and  p_ known,  $(^),  requires  only 
the  calculation  of: 

Z^  =  (ACL-np)//np(l-p) 

and  $(Z^)  =  Y  ^  cumulative  (left-hand)  normal  curve  table. 


(7) 


Arcsln  Transformation 


The  arcsin  transformation 

y  =  sin’^  /(r  +  3/8) /(n  +  3/4) 


produces  a  random  variable,  which  is  approximately  normally  distributed. 
(See  Johnson  and  Kotz,  1969.)  Thus  a  normalized  random  variable,  Z,  produces 
a  statistic,  the  asymptotic  distribution  of  which  is  normal,  where: 

=  2/ir[sin‘y^-^-|^  -  sin"^  /p]  (8) 

and  $(Z  )  =  Y  is  then  found  on  a  cumulative  normal  curve  table. 

Y 


PROBLEM  SOLUTION  -  EXACT  BINOMIAL 


Restating  the  problem,  the  objective  is  to  find  an  (j^,£)  pair  such  that: 


minimize:  n,c 

subject  to: 

P[r  <  c|n,pj]  > 

1-a 

(3) 

P[r  <  c|n,Pj,]  < 

P 

(4) 

where: 

p^  <  p^  and  p  < 

1-a. 

(5) 

Guenther  (1969)  develops  a  search  procedure  for  finding  (n_,£)  pairs  that 
satisfy  equations  3,  4,  and  5*.  While  Guenther's  algorithm  is  aimed  at 


♦Actually  Guenther's  paper  is  devoted  to  finding  sample  sizes  (n)  and 
acceptance  numbers  (c)  for  single  sampling  acceptance  plans.  The  procedure, 
however,  is  the  same.  In  his  paper,  the  hypergeometric,  binomial,  and  Poisson 
distributions  are  used. 


deriving  plans  by  hand  calculation  and  the  use  of  tables,  it  is  an  iterative, 
brute-force  technique  more  amenable  to  computerization  than  to  hand 
calculation. 

Hailey  (1980)  programmed  Guenther's  algorithm  to  find  the  minimum  single 
sampling  acceptance  plan  satisfying  equations  (3)  and  (4).  His  paper  contains 
the  FORTRAN  IV  computer  code  for  deriving  plans  using  either  the  binomial 
distribution  or  the  Poisson. 

The  algorithm  operates  basically  as  follows.  For  any  stipulated  value  of 

£,  there  is  a  minimal  sample  size,  satisfying  equation  (4).  That  value  is 

designated  a  minimal  value  of  For  the  same  value  of  there  also 

exists  a  maximum  value  of  satisfying  equation  (3).  That  value  is  designated 

n  ,  a  maximum  value  for  n.  If  the  solved  value  of  n  .is  less  than  the  solved 
—ft  —  —ft 

value  of  n  for  fixed  c,  no  feasible  solution  exists  for  that  value  of  c  or 
.ny  lesssf  value.  If  is  greeter  then  (or  equal  to)  n,.  then  any  value  of  n 
(a  unique  solution)  is  the  range  ^  <11  <  J1  feasible  for  that  value  of 
Cj,  In  fact,  feasible  solutions  exist  for  any  value  of  ^  greater  than  the 
designated  value,  as  well.  That  is,  an  infinite  number  of  plans  exist 
satisfying  equations  (3)  and  (4). 

The  search  procedure  begins  by  setting  equal  to  zero  and  solving 
for  n  and  n^  .  The  value  of  c_ is  increased  by  one  and  the  process  repeated 
until  a  feasible  range  of^is  found.  Hailey's  program  immediately  selects 
the  minimum  ji_,  jri^,  and  terminates  with  a  series  of  output  options.  A 
variation  of  this  program  was  used  to  derive  sample  sizes  and  ACL's  based  on 
the  binomial  distribution.  The  sample  size,  was  set  equal  to  ^  and  the 
Acceptance  Control  Limit 


ACL  =  c  +  0.5 


in  order  to  avoid  any  confusion  in  the  intepretation  of  points  falling  on  the 
control  limit. 


PROBLEM  SOLUTION:  NORMAL  APPROXIMATIONS 

Hand  calculation  using  the  algorithm  stated  in  the  previous  section  would 
become  most  tedious  and  time-consuming.  Available  tables  of  the  binomial  may 
not  cover  the  ranges  of  ^  or  required.  To  evaluate  the  binomial  where  c^ 
equal  50  requires  the  calculation  and  summing  of  51  terms,  which  is  a  large 
task  even  with  the  aid  of  a  sophisticated  pocket  calculator  or  small  com¬ 
puter.  This  procedure  would  have  to  be  repeated  many  times  before  the  minimum 
(n^,£)  pair  are  found.  However,  approximations  usually  require  the  evaluation 
of  only  two  equations  one  for  n_ and  one  f or 

Standard  Normal  Approximation 

If  denotes  the  value  that  cuts  off  an  upper  tail  area  of  o  under  the 
standard  normal  curve,  as  illustrated  in  Figure  2,  then  the  acceptance  plan 
(^,£)  pair  can  be  found  from  the  following  equations: 

c  -  n  pj 
/n  pj  (1-Pj) 

c  -  n  p 

-Z,  „  =  . -  ^  — 

/n  pg  (l-Pg) 


Solving  these  equations  simultaneously  yields: 


Figure  2.  Acceptance  Control  Charting 
Scheme  for  Binomial  Counts. 


The  values  of  and  may  be  calculated  rather  quickly  by  hand  or  by  the  use 
of  a  programmable  pocket  calculator.  This  is  the  simplest  formulation  of  the 
problem.  The  value  of  of  course,  is  rounded  to  the  nearest  integer. 

Arcsin  Transformation 

The  acceptance  control  chart  plan  (n.,£)  pair  can  be  found  by  solving  the 
following  equations; 

Zj.^  .  2 /S  [sla-l/HH  - 

-Zj.p  .  2  /F  -  '''>■*  -^23 

Solving  these  tvra  equations  simultaneously  yields: 

"  *  -sin"'  /P])3|^  (11) 

c  •  (n  ♦  3/4)  {sin  /")  ♦  sin"'  -  3/8  (12) 

Again,  the  value  of  jj^is  rounded  to  the  nearest  integer.  The  derived  function 
sin"^  (x)  «  tan'^  (x/ /l  -  )  may  be  used  to  obtain  the  inverse  sin  of  x, 

where  -1  <  x  <  1,  in  computers  which  do  not  have  an  inverse  sin  function, 
(Standard  Mathematical  Tables,  1973). 


ANALYSIS 


The  three  methods  of  calculating  (_n_,£)  pairs  for  Acceptance  Control 
Charts  were  investigated  using  the  popular  risk  levels  of  a  equals  0.05 
and  p  equals  0.10.  The  value  of  »*anges  from  0.005  to  0.06  in  increments  of 
0.005.  The  value  of  ^  changes  iteratively  in  order  to  maintain  specified 
levels  of  the  discrimination  ratio,  D_=22/£i.  The  value  of  Deranges  from  1.5 
to  5.0  in  increments  of  0.5. 

Acceptance  Control  Chart  plans  in  terms  of  (n_,  £)  pairs  are  presented  in 
Table  1.  The  three  methods  are  labeled  BINOMIAL,  NORMAL,  and  ARCSIN  to  iden¬ 
tify,  respectively  the  cumulative  binomial  distribution,  standard  normal 
approximation,  and  normalized  arcsin  transformation. 

Subgroup  sizes,  the^L's.  obtained  from  the  two  approximations,  ARCSIN 
and  NORMAL,  have  been  rounded  to  the  nearest  integer.  The  acceptance  control 
limits,  then's,  where  not  converted  to  integers  because  values  from  these 
approximations  yield  different  risk  protections  depending  on  whether  a 
continuity  correction  is  added  or  subtracted  from  before  it  is  converted. 

For  example: 


Acceptance 

Control  Limit  Truncate 

(Integerized)  Value  of 


Risk 

Protection 

Favored 


c'  =  c  +  0.5 

c  =  c 

c*  =  c  -  0.5 


Type  I  error 
Inconsistent 
Type  II  error 


where  0.5  is  a  continuity  correction  implicit  in  Laplace's  Theorem.  (See 
Johnson  and  Kotz,  1969,  p.53.)  Thus,  the  conversion  of  c_ to  integer  is 
performed  after  determining  the  continuity  correction  which  will  yield  the 
type  of  risk  protection  most  favorable  to  the  quality  control  program. 
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Tabl#  1  (b) 


PI"  0.0350  0.0400  0.0450 


Tabla  1  (e> 


Tablff  1  (d) 


L'  ■■ 
1.- ' 


I 


For  all  three  methods,  and  ^  decrease  as  the  difference  between  the  AQL 

(p^)  and  RQL  (£2)  increases.  This  can  be  observed  by  either  maintaining  a 
constant  discrimination  ratio  (D^  =  P2/P1)  and  increasing  pj,  or  by  maintaining 
a  constant  and  increasing 

Tables  2,  3,  and  4  indicate  the  risk  over-protection,  or  under-protection 
(negative  values),  which  each  plan  provides.  The  table  values  are  the 
differences  between  actual  (PCk<c|p])  values  and  stipulated  (l-o  or  p)  risk 
protection  probabilities.  The  actual  probabilities  were  obtained  using  a 
FORTRAN  program,  in  single  precision  mode,  on  a  POP  11/34  computer.  The 
probabilities  obtained  from  this  program  compared  with  published  binomial 
probability  tables  [N.B.S.  (1950),  Romig  (1953),  and  Harvard  (1955)]  with  an 
0.00001  accuracy. 

Binomial  Oistribution 

As  evidenced  in  Table  2,  ACC  plans  obtained  from  the  binomial  distribu¬ 
tion  over-protect  both  the  producer  and  consumer  by  having  low  probabilities 
of  committing  Type  I  or  Type  II  errors;  i.e. ,  rejecting  good  quality  products 
or  accepting  bad  quality  products.  The  largest  difference  between  actual  and 
stipulated  probabilities  is  2.21%.  The  over-protection  is  due  to  the  fact 
that  since  n^ and  c_ are  both  integers,  an  OC  curve  usually  cannot  pass  exactly 
through  points  l-o)  and  (2^,  p),  thus  the  search  procedure  finds  the 
minimum  2^ and  £  which  satisfy  inequalities  (3)  and  (4). 

Standard  Normal  Approximation 

The  acceptance  control  plans  obtained  from  the  standard  normal  approxima¬ 
tion  generally  have  smaller  subgroup  sizes,  n_,  than  the  plans  obtained  from 
the  other  two  models.  As  evidenced  in  Table  4,  this  method  consistently 
protects  against  Type  II  error  when  c*  is  used,  but  the  underprotection 
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Table  4.  ARCSINE  transforaation  risk  protection  difference  (R)  between  actual 
(P[k<c|p])  and  stipulated  (1-a  or  B)  probabilities  of  acceptance. 
Plans  In  Table  1. 


against  Type  I  error  may  be  as  large  as  14%,  When  ^  is  adopted,  protection 
against  Type  I  error  is  usually  attained,  or  is  at  most  1.25%  below  the 
stipulated  producer's  risk  probability,  but  under-protection  against  Type  II 
error  may  be  as  large  as  10%.  When  no  continuity  correction  is  used  with  c_, 
there  is  no  consistent  risk  protection,  and  under-protection  against  either 
type  of  error  may  be  as  large  as  6%. 

Arcsin  Transformation 

As  evidenced  in  Table  5,  acceptance  control  plans  obtained  from  the 
arcsin  transformation  consistently  provide  one-tail  protection.  When  c*  is 
used,  protection  against  Type  II  error  Is  attained,  but  not  against  Type  I 
error.  Under-protection  against  Type  I  error  may  be  as  large  as  9%.  When  c^ 
is  adopted,  protection  against  Type  I  error  is  attained,  but  under-protection 
against  Type  II  error  may  be  as  large  as  10%.  When  no  continuity  correction 
is  used  with  no  consistent  protection  is  attained,  and  under-protection 
against  either  type  of  error  may  be  as  large  as  3%. 

SIMULATION  STUDY 

A  computer  program  simulating  item  manufacture  and  control  charts, 
written  by  Davis  (1977),  was  used  to  analyze  the  performance  of  ACC's  derived 
from  the  binomial  distribution,  standard  normal  approximation,  and  arcsin 
transformation.  Twenty  replications  of  this  simulated  process  were  made, 
using  common  random  numbers  for  variance  reduction. 

Process  Description 

It  is  desired  to  have  producer  and  consumer  risks  of  5%  and  1D%, 


respectively.  The  cost  of  a  Type  I  error  is  considered  to  be  greater  than 
that  of  a  Type  II  error  in  this  simulation;  thus,  is  used  with  the 


Table  5.  BINOMIAL  ACC  Plan  Simulation  Results 
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approximations  because  of  Its  Increased  Type  I  error  protection.  In  Table  1, 
when  an  APL  of  1.5%  and  a  discrimination  ratio,  of  3.5  are  chosen  (I.e. , 
RPL  »  5.25%),  the  acceptance  control  plans  are  as  follows: 


Method 

n 

c 

c^ 

c* 

Plan  Chosen 

BINOMIAL 

175.0 

5.00 

(175,  5) 

ARCSIN 

182.0 

6.73 

7 

6 

(182,  7) 

NORMAL 

168.0 

5.11 

5 

4 

(168,  5) 

Thus  the  Acceptance  Control  Limits  were  positioned  at  5.5  for  the 
binomial  and  standard  normal  approximation  and  at  7.5  for  the  normalized 
arcsin  transformation.  The  simulation  "manufactured"  10,500  Items.  Initially 
the  manufacturing  process  had  a  nonconforming  percentage  of  1.5,  the  AQL. 

After  the  8,750th  Item  produced,  the  process  shifted  to  5.25%  nonconforming, 
the  RQL  level.  Thus  there  was  a  shift  from  the  expected  acceptable  process 
level  to  the  rejectable  process  level  to  test  the  plans  at  the  stipulated 
extremes. 

Control  Chart  Analysis 

The  control  charts  In  Figures  3  and  4  show  the  number  of  nonconforming 
units  found  In  each  subgroup  of  the  three  acceptance  plans.  Tables  5,  6,  and 
7  contain  the  simulated  data  used  to  plot  these  Figures.  The  process  shifted 
to  the  RPL  after  the  50th  subgroup  for  the  exact  binomial  plan,  and  In  the 
middle  of  the  49th  and  53rd  subgroups  for  the  arcsin  transformation  and  normal 
approximation  plans,  respectively. 

The  binomial  plan  made  one  Type  I  error  (subgroup  20)  In  50  subgroups  of 
good  quality  and  no  Type  II  errors  were  made  In  the  final  10  subgroups.  The 
arcsin  transformation  plan  made  one  Type  I  error  (subgroup  32)  out  of  48  good 
quality  subgroups  and  no  Type  II  errors  In  the  last  ten  subgroups.  The 
standard  normal  approximation  plan  made  no  Type  I  errors  In  the  first  52 
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Figure  3.  Acceptance  Control  Chart,  Exact  Binomial 


Process  y,  Process 

In-Control  1^,  Out-of -Control 


Table  6>  ARCSINE  ACC  Plan  Slnulatlon  Results 
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Table  7.  NORMAL  ACC  Plan  Simulation  Results. 
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subgroups,  but  had  three  Type  II  errors  (subgroups  57,  58  and  60)  in  the  last 
eleven  subgroups,  when  the  process  had  shifted  to  the  RPL. 

The  total  number  rejected  in  a  subgroup  may  be  expected  to  exceed  the  ACL 
for  either  of  two  reasons,  (1)  the  existence  of  assignable  causes,  or  (2)  the 
existence  of  a  quality  level  which  exceeds  the  APL.  In  either  case,  the  only 
clue  given  by  the  acceptance  control  chart  as  to  the  cause  of  lack  of  control 
Is  the  time  at  which  lack  of  control  at  the  desired  level  was  observed.  For 
this  reason.  Immediate  corrective  action  should  be  taken  whenever  a  point 
exceeds  the  ACL.  This  simulation  did  not  Include  such  corrective  action 
because  It  was  desired  to  observe  the  consistency  of  each  control  chart  In 
providing  the  desired  risk  protections. 

CONCLUSIONS 

Three  methods  for  obtaining  subgroup  sizes  and  acceptance  control  limits 
were  compared.  In  addition  to  utilizing  the  exact  binomial  distribution,  the 
standard  normal  approximation  and  a  normalized  arcsin  transformation  were 
used.  Acceptance  plans  obtained  by  using  the  binomial  distribution  provide 
the  stipulated  risk  protections  (guaranteed  over  protection)  for  both  producer 
and  consumer.  This  results  from  the  strict  application  of  the  Inequality 
constraints  of  equations  (3)  and  (4).  When  a  continuity  correction  of  0.5  Is 
added  to  (subtracted  from)  the  ^derived  from  either  of  the  latter  two 
approximation  methods,  strict  risk  protection  against  Type  I  (Type  II)  error 
Is  attained,  but  not  against  both  error  types.  If  no  continuity  correction  Is 
used,  the  risk  protection  from  these  approximations  Is  Inconsistent;  I.e. , 
risk  protection  alternates  against  both  types  of  errors,  with  no  discernable 
pattern. 
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For  the  risk  levels  studied  (1-a  >  0.95  and  p  «  0.10)  and  the  wide  range 
of  values  of  21  and  ^studied,  the  normalized  arcsin  transformation  yielded 
results  closer  to  design  than  did  the  standard  normal  approximation.  Risk 
protection  losses  ranged  as  high  as  2.96%  for  the  producer  (1-a)  and  3.38%  for 
the  consumer  (p)  with  no  continuity  correction  factor  applied.  These  are  the 
maximum  (underlined)  negative  numbers  In  Table  4.  With  a  continuity  correc¬ 
tion  factor  of  0.5  added  to  the  solving  value  of  c,  the  producer  received  at 
least  the  required  protection  but  the  consumer  loss  of  protection  reached  as 
high  as  9.81%  (nearly  doubled).  With  0.5  subtracted  from  the  solving  value  of 
c_,  consumer  protection  at  the  specified  level  was  assured  but  loss  of  producer 
protection  Increased  to  9.08%,  I.e.,  from  a  design  level  of  0.05  to  as  high  as 
0.1408.  Unless  protection  at  one  level  Is  vital,  as  opposed  to  protection  at 
the  other  level,  no  continuity  correction  Is  recommended  when  the  arcsin 
transformation  Is  to  be  used. 

Results  from  the  standard  normal  approximation  were  not  as  good.  In 
general,  as  those  achieved  by  applying  the  arcsin  transformation.  Without 
adjustment  by  a  continuity  correction  factor,  loss  of  producer  protection 
ranged  as  high  as  6.92%  and  loss  of  consumer  protection  as  high  as  4.29% 
(negative  underlined  values  In  Table  3).  With  0.5  added  to  the  solving  value 
of  Cf  loss  of  producer  protection  was  reduced  to  1.2%  but  loss  of  consumer 
protection  was  Increased  markedly  to  10.19%.  With  0.5  subtracted  from  the 
solving  value  of  c,  loss  of  producer  protection  Increased  to  13.64%  but  con¬ 
sumer  protection  at  the  design  level  was  assured.  As  was  the  case  with  the 
normalized  arcsin  transformation,  no  continuity  correction  can  be  recommended 
unless  It  Is  Imperative  to  meet  (or  nearly  meet)  the  design  level  of  protec¬ 
tion  at  either  the  producer  or  the  consumer  quality  protection  level. 


Both  the  standard  normal  and  arcsin  transformation  require  the  evaluation 
of  only  two  equations  to  obtain  an  ACC  plan  (n.*£)  pair.  These  equations 
easily  may  be  evaluated  using  Standard  Mathematical  Tables  or  a  scientific 
calculator,  since  they  only  require  the  use  of  square  root,  sine,  and  Inverse 
sine  functions.  The  binomial  distribution  requires  at  least  c_  evaluations  of 
equation  2,  which  contains  factorial  and  exponential  terms,  and  a  search  for 
the  minimum  n_  among  various  (n^,£)  pairs  that  satisfy  the  Inequalities  In 
equations  3  and  4.  However,  even  most  home  computers  have  the  capability  of 
performing  these  evaluation  and  search  tasks  quickly;  they  would  be  tedious 
and  time  consuming  If  performed  using  tables  and/or  pocket  calculators. 

Finally,  to  assure  the  desired  protection  for  both  producer  and  consumer, 
the  exact  binomial  should  be  used  to  obtain  ACC  plans,  provided  computer 
facilities  are  available.  The  normalized  arcsin  transformation  Is  preferable 
to  the  standard  normal  approximation  because  Its  likely  degree  of  under- 
protection  Is  about  half  that  of  the  standard  normal.  Possible  under-protec- 
tlon  afforded  by  the  standard  normal  Is  about  double  that  of  the  arcsin  In 
absolute  terms. 

A  complete  listing  of  the  computer  programs  used  to  develop  and  evaluate 
the  various  plans  Is  provided  In  the  Appendix. 
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^K  a^  ^a  ^U  ^U  ^U  ^U  ^U  ^U  ^b  ^U  ^U  «^  ^U  ^U  ^J  a^  a^ 

^\  1^  I '  1  ^p  ^p  ^p  ^p  ^p  ^p  ^p 

20  REM*  program:  BINOMIAL  ACCEPTANCE  CONTROL  PLAN  (N»C)  * 

w\^  rvir. 

'^O  niM  SUMLOG(IOOO) 

50  NMAX  ^  1 
60  SUMLOG<NMAX)  =  0 

70  INPUT  ’ENTER  PROHUCER  %  CONSUMER'S  RISKS  ’»  A»B 

SO  INPUT  ’ENTER  ACCEPTABI.  E  8  REJECTABI.  E  PROCESS  LEVELS  ’»  PlfP? 

90  OOSUP  170 

100  PRINT  ’SUBGROUP  SIZE  N  ’  ;N 
110  PRINT  ’CONTROL  LIMIT  C  =  ’JC 
120  LNIi 

1**^  /\  1^  1^*  1^  ^U  a^  ^a  ^U  a^  a^  a^  ^a  a^  a^ 

•  \r  fT«  I  1  ^P  ^p  ^P  ^P  ^P  ^P  ^P  ^p  ^P  ^P  ^P  ^P  ^P  ^P  ^P  ^P  ^p  ^P  ^P  ^P 


140 

REM* 

subroutine: 

SINGLE  SAMPLING  PLAN  SEARCH 

* 

1 50 

REM* 

* 

160 

REM* 

REFERENCE : 

GUENTHER  (1969)  8  HAILEY  (1980) 

* 

1 70 

REM* 

* 

1 80 

REM* 

given: 

A  == 

PRODUCER'S  RISK 

* 

1 90 

REM* 

B  = 

CONSUMER'S  RISK 

* 

200 

REM* 

PI  = 

ACCEPTABI.E  QUALITY  LEVEL  (AQL) 

* 

210 

REM* 

P2  = 

REJECTABI  E  QUALITY  LEVEL  (RQI  ) 

* 

220 

REM* 

*  * 

230 

REM* 

finds: 

N  = 

MINIMUM  SUBGROUP  SIZE 

* 

240 

REM* 

C  = 

NONCONFORMANCE  CONTROL  LIMIT 

* 

250 

l"a  IT'  4^ 

I*  J  ^4  ^P  ^P  ^p  ^P  ^P  ^P  ^P  ^p  ^p  ^P  ^p  ^p  ^p  ^p  ^P  ^P  ^P  ^P  ^P  ^P  ^P  ^P  ^p  ^p  ^p  ^P  ^P  ^p  ^P  ^p  ^P  ^P  ^P  ^p  ^p  ^p 

260 

N  = 

1 

270 

c  =  - 

1 

280 

C  =  C 

+1 

290 

P  =  P 

<5 

300 

N  =  N 

"+1 

310 

GOSUB 

370 

\  REM***  CALL  BINOMIAL  (P2»  Nf 

C»  PACC) 

320 

IF 

PACC  >  B 

THEN 

300 

330 

P  = 

PI 

340  GOSUB  370  \  REM*#*  CALL  BINOMIAL  (P3f  N»  Cf  PACO 

350  IF  PACC  <  (1-A>  THEN  2S0 
360  RETURN 


370  REM***************************************************************** 


380  REM*  subroutine:  CUMULATIVE  BINOMIAL  PROBABILITY  * 
390  REM*  * 
400  REM*  given:  N  =  SUBGROUP  SIZE  * 
410  REM*  C  =  NONCONFORMANCE  CONTROL  LIMIT  * 
420  REM*  P  =  PROBABILITY  OE  NONCONFORMANCE  * 

430  REM*  * 
440  REM*  EINnS:  PACC  =  PROBABILITY  OF  ACCEPTANCE  (CUM.  BINOMIAI  )  * 


450  REM***************************************#********#**************»!t: 

460  0  =^  1-P 

470  CUMULA  •=  Q''N 

480  IF  C  ==  0  THEN  600 

490  IE  N  <>  NMAX  THEN  540  \  REM***  LOG  SUMS  AIRFAHY  IN  MEMORY 

500  FOR  K  =  <NMAX+1)  TO  N  \  REM***  COMPUTE  ONLY  NEW  LOG  SUMS 

510  SUML06<K)  =  L0G10<K)  +SUMI0G(K-1> 

520  NEXT  K 

530  NMAX  =  N  \  REM***  LARGEST  LOG  SUM  (FACTORIAL)  IN  MEMORY 

540  PLOG  =  LOGIO(P) 

550  QL.OG  =  I  OGIO(O) 

560  FOR  K  =  1  TO  C  \  REM***  COMPUTE  CUMJILATIVE  PROBABILITY 

570  FACTOR  =  SUMI.OG(N)  -SUMLOG(N-K)  ~SUMLOG(K) 

580  CUMULA  =  10" (FACT OR  +K*PI  OG  +(N~K)*OI  OG)  +CUMULA 

590  NEXT  K 

600  PACC  =  CUMULA 

610  RETURN 


S I NRI.  F;  » S I NPLE/L  in  =^S  I NGI.  H 

nr - 

nc  program:  compares  single  sampling  pi.ans  for  a  process 

CC  WHICH  HAS  BINOMIAL l.Y  BIS TRIBIJTEP  DEFECTIVES* 

CC 

CC  programmer:  CARLOS  AMAPO  fall  'SI 

CC 

CC  RESEARCH  FOR:  PR.  R*l.  *  LEAVENWORTH^  ISF  PEPT»  UNIV.  OF  FLORIDA 

CC - 

0001  DIMENSION  CLlJ('.?)r  CLL(?) 

0002  CALL.  ASSIGN  (2» 'SINGLF.DAT' >  !  OUTPUT  EII.E 

0003  130  CALI.  PUTCHR  <2A»S)  !  CLEAR  CRT  SCREEN 

0004  WRITE  <5fl31> 

OOOS  J31  FORMAT  ('OENTER  DATA  (START  IN  1ST  I.ETTER  OF  TITL.E)'// 

I  '  PO  ALFA  PI  BETA  ITERS') 

OOOA  READ  <f)»133TFRR=^130)  P0»  AI.FAf  Pit  BETA»  ETERS 

0007  133  FORMAT  ( (F7 . 4  ,  IX  )  ) 

0008  IF  (PO.LF.O)  GO  TO  199 

0010  IF  <P0.FQ*P1)  GO  TO  127 

0012  IF  (ETERS. NF*0)  GO  TO  124 

0014  GO  TO  129 

C 

OOlS  127  WRITE  (S»*) 

0014  WRITE  (f.f*)  'PO  CAN  NOT  E(?UAL  PI  (7ER0  DIVIDE  ERROR)' 

0017  GO  TO  194 

C 

0018  124  CALL.  ITERS  (ALFA*  P0»  BETA»  Pl»  Af  FTER8) 

0019  IF  (A.FQ.'Y')  GO  TO  194 

0021  IF  (A.EO.'R')  GO  TO  194 

0023  GO  TO  130 

C 

0024  129  WRITE  (2fl2S)  POf  ALFAr  Pl»  BETA 

002L=i  128  FORMAT  (///'  GIVEN  (P0»  Af  Pl»  B) : '  »4(F7.4r2X) ) 

0024  CALL  ARCSIN  (AI.FA^  P0»  BETA»  Pl»  SN»  Cr  CLU»  CLL) 

0027  WRITE  (St#)  SNt  Ct  '  FROM  ARCSINE' 

0028  WRITE  (2.T#)  '  ' 

0029  WRITE  (2r#>  'ARCSINE  TRANSFORMATION' 

0030  M='  1 

0031  GO  TO  134 

C 

0032  132  CALL.  NORMBI  (AI.FAr  POt  BETAt  PI  t  SNt  Ct  CLUt  CLL.) 

0033  WRITE  (2t*)  '  ' 

0034  WRTTF  (2t*)  'NORMAL  APPROXIMATION' 

003S  M=^  2 

0034  GO  TO  134 

C 

0037  134  CALL  GUNTHR  (AI.FAt  POt  BFTAt  PIt  SNt  C) 

0038  WRITE  (2f*)  '  ' 

0039  WRITE  (2t*)  'EXACT  BINOMIAL.' 

0040  M-^  3 

C 

0041  134  IF  (C.En.999.)  GO  TO  199 

0043  CALL  EXACT  ( SN t PO t C t BXL FC t N ) 

0044  CALL  EXACT  (SNtPI tCtPtN) 

C 


a3 


I 


siNni.Ef  siN(5i.e:/i..i :  irrsiNOi..F 


004  s 
004A 

l.?55 

URITF  <2t135)  SN»  BXI.FCr  Hr  P 

FORMAT  ('  RIMOl.F  SAMPLINO  PI  ANJV/ 

0047 

0048 

194 

I  '  8AMPI.F  8I7F  ^'»F9.2»5X» 

I  '  P  >  1-A  =-'rF9.4/ 

I  '  MAX  nFFFCTS  ^'»F9.2»5X» 

I  '  P  <  B  =r'rF9.4) 

no  TO  (i5:?.i;?A»i30)»  m 

WRJTF  (5r#> 

0049 

OOSO 

1.9A 

WRITF  (Sr#)'  INPUT  FRROR ‘  TYPF  “R*  TO  RF8TART 
RFAB  (5rl97)  A 

0051 

197 

FORMAT  (Al) 

0052 

0054 

1.99 

IF  (A.FO. 'R' .0R»A.F(J4 'NM  GO  TO  li^O 

WRITF  <5f*) 

0055 

005A 

STOP  '  HAOF  A  ROOD  1.  IFF' 

F.NJ) 

a4 


BINRI.  F. » SINRI..F/I..1 : 1  =SINfil.  E 


0001 


StJBROiniNI:  RIJNTHR  (ALFAr  P0»  BETAf  Pl»  SN»  C> 


prorram:  obtains  exact  binomial  sinri.e  sampi.inr  plan 

RFFt  01IE:NTHFR»  W.C*»“IJSE  OF  BINOMTALr  HYPFRRF OMI;  TRIO r  ANB 

POISSON  TABI  ES  TO  OBTAIN  SAMPI.INR  PLANS'^  JOURNAL  OF 
OIJAI.  ITY  TE:CHt  UOI.  Ir  NO.  2f  APRIL  IVAVr  PP  lOS-109. 


n=  -1.  !  START  AT  OO 
SNS:-^  1,  !  START  AT  SNS=2 

C=  r+  1  . 


SNS=  SNS+  1. 

CAl.l  EXACT  <SNS»  Pir  Cf  BXI.EC»  N) 

IE  (bxlec.gt.be:ta)  ro  to 

SNI  SNS-  1.  !  SMALLER  SNL  NOT  NEEBFn 
SNL-  SNI.+  1. 


WRITE  <fjr*)  SNSr  C»  SNI.  »  BXLEC 
CALL  EXACT  <SNI.  »  P0»  V.s  BXI.ECf  N) 

IE  (BXLEC.Rr.Cl .-ALFA))  RO  TO  4 

SNI.=  SNI  -  1.  !  PREUIOUS  SNI..  IS  THE  ONE  WANTED 

IF  (SNS.RT.SNL)  RO  TO  2 

SN"=  SNS 

RETURN 


SINBI. 

0001 


ooo:? 

0003 

0004 


0005 

OOOA 


0008 

0010 

0011 


0012 

0014 

0015 
001 A 
0017 
0018 


0019 

0020 


FfSINRI.F/l.  mr-SINGLF 


nc 

nc 

cc 

nc 

nc 

rc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc 

cc- 

cc 

cc 

cc 

cc 

cc 

cc 

cc- 


c 

c 


c 

c 

c 


100 

c 

c 

c 


110 


1 1 1 
c 
c 
c 


c 

211 


c 


SUBROUTINF  FXOCT  (SNt  P»  C»  BXI.FCr  N> 


program:  COMPUTFS  PXACT  comm,  binomial  PROBABIl.  ity  of  x.i.f  ,c 


mf:thod: 

A4  “N"  number  of  I  00(10)  SUMS  ARF  COMPUTFH  ONCF,  ONI.  Y » 
AND  STORFB  IN  “SUMI..On(  I  >  •  OFCTOR. 

B.  10**(  BUMI.OOd)  >  RFPRFSFNTS  I-FACTORIAI.  1  THFRFFORF» 
FACTORIAL  MUI. TIPI.  ICATIONS  ARF  RFDUCFTi  TO  8UMKATI0N8» 
Al  l.  FACTORIALS  NF.FBFB  ARF  IN  STORAOF  UP  TO  SUMIOG(N), 

C.  WHFNFVFR  ITFRATIONS  ARF  RUN»  THF.  SUMI..00 ( I ) ' 8  IN  MFMOI 
ARF  NOT  RFCOMPUTFIi . 

n.  N  >  C  >  0?  N-FAC  >  C-FAC;  ANB  SUMI..OO(N)  >  8UMI  00<C). 
F.  SUMl.OOd)  SUBTRACTION  ==•  I-FACTORIAL  rUOISlON. 

F«  1.00(10)  Al.l.OWS  FFFICIFNT  HANIH  INO  OF  I.AROF  NUMBFRS. 


PROGRAMMFR:  CARI.OS  H.  AMAHOf  ISF;  nFPT.f  UNIO.  OF  FI.ORiriA 


OIVFN:  SN  SAMPLF  NUMBER  (SIZF) 


P  -  PROB.  OF  DFFFCTIVFS 

C  =  ♦  OF  BFFFCTIVFS  IN  SAMPLE 


computes:  BXLFC  binomial  PROB  (X.I  E.C) 

SUMLOGd)  ==  VECTOR  CONTAINING  SUM  OF  LOO(l)  THRU  LOO 


VIRTUAL,  SUMI.  00(9000) 
0  l.-P  ' 

NN  -  8N 


»>  BINOMIAL  PROB,  WHEN  C^^O  «< 

CSUMS  =  0  *#NN 

IF  (C^EQ.O.)  00  TO  333 

»>  AVOID  RFCOMPUTINO  SUMLOO(I)'S  ALREADY  IN  MEMORY  «< 

IF  (NN. or, 9000)  00  TO  998 
IF  (N-NN)  100»211f2il 
M:--^  NT  1 

»>  COMPUTE  N  SUMLOOS  — I:  OUI VALL  FNT  TO  N-FACTORIAL.  <« 

IF  (M.OT.l)  00  TO  110 
SUMl.OOd)  ^  0, 

M  2 

DO  111  1=^  M»  NN 

SUMl.OOd)  ^  Al  OOlO(FLOATd)  )  +SUMLOO(T-1) 
CONTINUE 


»>  COMPUTE  C  CSUMS  — FOUIVAI.LFNT  TO  SUM  OF  PROB  COMBINATIONS  4 
T.F. f  CUMMULATIVF  BINOMIAL  DISTRIBUTION  COMPUTATION 


KC  -  C 
N^^  NN 


a6 


RINRI  H  fSINOLF/I  I!  ]=-BINOI.F 


r; 

c 

»>  SUM  (ftCCUMULATE)  PROBABII.  ITIF  8  <« 

0021. 

PI.OR=-  AI..0G10<P) 

0022 

01.08--  Al. 081  0(0) 

002;? 

no  322  K=--  1»  KC 

002/1 

CSUMS  =  1.0. *#(  SUMI..08<N)  -SJJMI. 08<N-K )  -SUMI.08(K) 
Z  +K*PL00  +<N-K)#0I.08  )  +OSUMS 

002rj 

322 

CONTINUrt 

0026 

U 

333 

CSUMS 

0027 

RETURN 

0028 

998 

WRITF  (Hr 999)  SNr  Pr  C 

0029 

999 

FORMAT  ('OWHOOPS!!!  SAMPl.F  SIZE  -'rFlO.2// 

7.  ‘  AT  P  -'»FT.0.4»'  AND  C  -"r  F1.0.4) 

n  n  rsn  n 


s  I  NBi.  F » s  I  Noi.  f;/i.  i;i=^single 


0001 


ooo:? 

ooo;? 

OOO'l 

ooos 

0006 


0010 

0011. 

ooi;t 

0014 

001  s 

OOIA 

0017 

OOIB 

0019 

oo:?o 

0021 

0022 

0022: 

0024 


002A 


SUBRnUTINI;  NORMBI  <AI  FA»  P0»  BFTAt  Plf  SN»  C»  CLUf  CIJ.  ) 

- - 

c  program:  computes  a  sinri.e  sampi.img  pi.an  using  the 

C  NORMAl  APPROXIMATION  FOR  A  BINOMIAL  BISTRIBUTION 

C 

C - 

DIMENSION  CI.U<2)i  Cl.l  (2)»  ZV(12>»  7T<12)»  7(2) 

DATA  7V(l)/2.;?2A/»  7V<2>/2.054/f  7V(7)/1 .8Bl/»  70  (  4  ) /I  ♦  7S1/ » 
I  7V(5)/1  .A45/»  7V(A>/1  .SSS/f  70(  7 ) /I  .  474/ »  70  ( B ) /1 . 40fj/ » 

I  7V<9)/1  .;’f4/f  7V<10>/1  .2B2/»  70 ( 1 1  ) /1 . 02A/ »  70(12)/. B42/ 

DO  29  1»  10 

29  ZT(I)=  1.-  I/IOO. 

ZT(11)=  .BS 

ZT(12)=  .BO 
P7^  l.-ALFA 
1 

30  DO  31  1=^=  It  12 

IF  (P7.NF.7T(I))  GO  TO  31 
7(K)=  70(1) 

GO  TO  (32t33)t  K 

31  CONTINUn 

WRITE  (5t*)  'error:  ALFA  OR  BETA  NOT  IN  DATA  STATEMENT' 

C"  999. 

RETURN 


32  PZ:^  l.-BETA 
K~  2 

GO  TO  30 


33  POO"  (1-  P0>*  PO 

pq)=  (1-  pi)»  pi 

SN-  ((Z(l)»  SORT(PaO)  +  7(2>*  SQRT(Pni))/  (PI-  P0))*«2 
C=  7(1)*  SORT(SN*PnO)+  SN*  PO 


RETURN 


SINBlF  f  SINBI..b:/I..I ;  1=^=SINGI.F 


0001  SUBRruniNir  arcs  in  (AI.FAt  P0»  BETA»  Rl»  SNr  C) 

CC - 

CC  PROPRAM:  SINGI.F  SAMPLINB  plan  using  an  ARCSINir  NORMALIZING 
nr  TRANSFORMATION*  TO  APPROXIMATF  A  BINOMIAL  PR0CF8S. 

rc 

CC  1.  NORMAL  ARC8TNI:  TRANSFORMATION 

CC  2,  FRFFMAN  X  TUCKFY  (F  X  T) 

CC 

CC  RFF:  JOHNSONf  Nf  S.  KOTZr'BIST.  IN  STATISTICS  -niSCRFTH  DIST.'i 
CC  HOUGHTON  MIFFl.IN  COf  BOSTON t  19A9r  P  AS. 

CC 


CC 

giofn: 

ALFA 

PROBUCFR'S  RISK  PROBABILITY  OF  REJECTION 

CC 

BFTA  ^ 

CONSUMER'S  •  •  •  ACCEPTANCE 

CC 

PO 

ACCEPTABLE  PROCESS  LEVEL 

CC 

PI 

REJFCTABLE 

CC 

CC 

computes: 

SN 

SAMPLE  SIZE 

CC 

C 

MAX  NUMBER  OF  BErFECTIVES  IN  ACCEPTANCE 

CC 

CC 

VARIABLES 

• 

4 

CC 

7(1) 

-  7(1~ALFA) 

CC 

7(2) 

=  7.(1 -BETA) 

CC 

SINV(l) 

ARCSINE  OF  SORT (1 -ALFA) 

CC 

SINV(2) 

-  -  •  ■  (1-BETA) 

CC 

SS(1) 

-  SAMPLE  SIZE  OF  NORMAL  ARCS  INI:  TRANSFORMAT  10 

CC 

SS(2) 

_  .  ■  •  F  X  T  ■ 

CC 

B(l) 

^  ACCEPTABLE  ♦  BEFECTIVES  IN  NORMAL 

CC 

CC— 

B(2) 

*  ACCEPTABLE  ♦  BEFECTIVES  IN  F  X  T 

0002  BIMITNSION  Z<2)»  SINM(2)»  SS(2)f  0(2) »  7T<12>»  ZV(12)»  P(2) 

0003  BATA  7V<1  )/2.32A/f  7.V<2.>/2.0S4/»  7V(3)/1 .881/f  7U(4>/1 .751/» 

I  7V<S)/l.A45/»  ZM(6)/l.SSS/»  7V(7>/1 .47^/»  7V<8>/1 . 40S/ » 

I  7V(9)/1.34/*  7M<10)/1.282/»  7V( 1 1 >/l .03A/»  7V(12)/.842/ 
0004  BO  29  I^  1»  10 

OOOS  29  7T<I)=  1,-  I/IOO. 

OOOA  7T<ll):r-  ,BS 

0007  7T(12)=  .80 

0008  P(l)"  PO 

0009  F<2}^  PI 

C 

0010  P7:^  1.  -  AI.FA 

0011  K-  1 

0012  30  no  31  I-  If  12 

0013  IF  (PZ.NF.ZTd))  GO  TO  31 

001  S  7<K)==-  70(1) 

001  A  SORP^  SORT(P(K)) 

0017  SINV(K)=  ATAN(  S«RP/  S«RT( -S«RP#SQRP+1 )  )  !  ARCSIN 

(>(>18  GO  TO  (32r33)f  K 

0019  31  CONTINUI: 

0020  WRITH.  (Sf»)  'FRRORJ  ALFA  OR  BFTA  NOT  IN  BATA  STATFMFNT' 

0021  (>-  999. 

0022  RHTURN 


S I NRI.  F. » S I  NBI..F./ 1..  in  NGI.F 


oo:?-i  K=  2 

002S  GO  TO  30 

r 

00?A  33  SB<?>=  (  (7(J.)+7<:?)>/  (SINV<2)-SINV(  1  )  )  )**2 

00:?7  SS<1)=  SR<?)  /A 

00:?B  X=  S1NV(J.)+  7(1)/  (:>#B0RT(.SS(1  )  )  > 

0079  0(1):^  (SSd  )+0.7f:)  ♦  (BIN(X>*#7)  -(3./B.) 

n 

0030  BN=  SS(l) 

0031.  i>-  rin> 

0037  RFTURN 


SINRI.  t  ..E/l.  I  i  1=SINRLF 


0001 


0002 

ooo;? 

0004 

ooos 

OOOA 

0007 

OOOB 

0010 

0012 

0013 

0014 

001 A 
0017 
0018 
0019 
0020 
0022 

0073 

0024 

002A 

0027 

0078 

0030 

0031 

0032 

0034 

003S 

003A 

0038 

0039 

0040 

0041 

0042 

0043 

0044 
004  Fi 
004  A 
0047 
0048 
0049 
OOSO 
0051 


SUBROiniNir  n>:R8  <AI.Fft»  P0»  BFTA»  P1»  Af  FTFR8) 

cc - - - 

nn  program:  bftfrminfs  bfvfrai  singlf  sampling  pi.ans 

nr  BY  nHANRlNG  PARAMFTFR  VALUFS  IN  FACH  ITFRATION. 

nn - 

niMFNSION  88(3) »  08(3) »  B0(3)f  Bl(3) 

COMMON  /FTFR/  Al.  » A8t  AHr BL»BS»BH»POI..f  POSrPOHtPlI.  »P1S»P1H 
WRITF  <5rl01) 

101  FORMAT  ('OFNTER  ITFRATION  (STFP)  ANB  HIGHEST  OALDFS'// 

7  '  (ALFA)  ALFA  (PO)  PO  <P1>  PI  (BETA)  BETA' 

READ  (5»103)  AS»  AH»  POS»  POH»  PlSr  P1H»  BSr  BH 

103  FORMAT  <8F7.4) 

r 

IF  (AS.I  T.O.OR.AH.LT.AIFA)  GO  TO  194 
IF  <BS.GT*0)  GO  TO  104 
BB:^  AS 
BHr=  AH 
C 

104  IF  (POS.I  T.O.OR.POH.l  T.PO. OR. PIS ♦  I  T. 0. 0R.P1H.lt. PI. OR. 

7  PIS. I  T.POS.OR.BH.I  T.BFTA)  GO  TO  194 

Al.-  ALFA 
BL=^  BETA 
POL^  PO 
Pll.=^  PI 

IF  <FTFR8.EQ.8. >  GO  TO  114 
GO  TO  no  !  BEGIN  ITERATIONS 
C 

lOA  Pl=  P1+  PIS 

IF  (Pl.LF.PlH)  60  TO  110 

Pl==-  PIL 
PO-  P0+  POS 

IF  (PO.LF.POH)  GO  TO  110 
P0=  POL 

ALFA-  ALFA+  AS 
IF  (AIFA.LE.AH)  60  TO  110 

ALFA-  Al. 

BETA-  BETA+  BS 
IF  (BETA.LE.BH)  GO  TO  110 
WRITE  <5f#) 

WRITE  <5»*)'  FNB  OF  ITERATIONS!  END  SESSION?  IIY/NI' 

A=  'Y' 

RETURN 
194  A^'R' 

RETURN 
D 

no  CALL  ARCS  IN  (ALFAr  P0»  BETA»  Pl»  8S(1)»  CC(1)) 

CALI  NORMBI  <ALFA»  P0»  BETAf  Pl»  8S<2)»  CCC?)) 

CALL  GUNTHR  (AI.FAf  P0»  BETAf  Plr  SS<3>»  CC(3)) 

DO  117  1=^  If  3 

CALL  EXACT  <SS(I)f  POf  CC(I>r  BO<I>f  N) 

117  CAl.l.  EXACT  (SS(I>f  Plr  CC(I)f  Bl<I>f  N) 

Al=^  l.-ALFA 

WRITE  <7»113)  Alf  BETAf  POf  Plf  SS<3)f  CC(3)f  P0(3)f  Bl(3)f 
7  SS<l)f  CC(l)f  B0(1>»  BKDf  SS(2)f  CC<7)f  B0(7)f  Bl(7) 


N  NN  Nj 


r>  I  Nfil.  K  »  S I NRI  E/L 1 : 1  =S  I  N(il.  F 


00S2  11.:?  FORMAT  ('OOIVEN  (1-A»  B»  P0»  PI  >  t  '  f  4 ( F7 . 4 »  '  f  '  ) 

'  SINOI.E  SAMPl  INO  Pl.ANt  N'»7XrM: 

'  EXACT  BINOMIAI.  '  f8X»2FB.2f2F8.4/ 

"  ARCSINE  TRANSFORMATION' f2FR. 2 »2F8. 4/ 

'  NORMAI..  APPROXIMATION  ' » 2F8 . 2 » 2F8 . 4  ) 
00  TO  10<^^ 

1.14  CALL  ITER88  (AI.FAr  POr  BKTAf  PI  r  A) 

RFTURN 


oos;? 

0054 

0055 


SINGI  F»f>TNGLE/LI :  l^^SINGI.E 


0001 


C==== 

c 

n=:==* 


0003 

0003 

0004 

0005 

113 

OOOA 

0007 

C 

0009 

0010 

114 

0011 

0013 

0013 

0014 

115 

0015 

OOIA 

C 

0017 

0018 

0019 

llA 

0020 

0033 

C 

0034 

0035 

003A 

C 

0027 

0029 

0030 

0031 

111 

C 

0032 

0033 

0034 

132 

0035 

003A 

1  17 

0037 

0038 

119 

0039 

0040 

131 

0041 

133 

0043 


SUBROUTINF  ITHRflfl  <AI.FA»  POr  BHTA»  Plf  A) 

routine:  compares JOFSCRIMINATION  ratios  (Pt/PO) 

DIMENSION  SS<9)t  nC<9),  B0(9>f  Bl<9)»  PACOr  PBC.?) 

COMMON  /ETER/  Al.  r  AS » AH » Bl.  » BS » BH  f  POI.. » POS  r  POH  f  P 1 1.  f  P 1 S  » P 1 H 
WRITE  (S»*) 

WRITE  (5»*) 'ENTER  D-STEP  FOR  PI' 

READ  (.fit*)  DS 
IF  <DS,l.T.O>  00  TO  11? 

1 

Al~  1,-AI.EA  ! 

WRITE  (:?fllS)  A1»AI..FA,PETA 
WRITE  (l»llf.)  Al» ALFA » BETA 

FORMAT  ('OOIVEN  (1-Ar  Ar  B > ! ' » 3 ( F7 . 4 » ' r ' ) > 

PA(1)=  POI.. 

PB<1)=  Pll.. 

PH=-  P1H+  0.001  I 

SP1=  PIS 

D=^  SPl/POS 
DO  111  M=  1»  3 

00==  PB(M)/PA(M> 

IF  (DD+O.OI.IE.D.OR.DD-O.OI.GE.D)  00  TO  138 
IE  (PB<M) .OT.PH)  00  TO  136 

CAl.l.  OUNTHR  <AI.EA»  PA(M>f  BETA»  PB(M>»  8S(M>»  CC(M)> 

CAI.I.  ARCSIN  <Al.FA»  PA(M>»  BETAr  PB(M)»  SS<M+3>»  CC(M+3> 

CALI.  NORMBI  <AI.FA»  PA(M)f  BFTA»  PB(M>»  SS(M4A>f  CC(M+<!*) 

IF  <M*0E.3)  00  TO  111 
PA<M+1)=  PA<M)+  POS 
PB<M+1)-  PB<M)+  SPl 
CONTINUE 

DO  133  M=’=  Xt  9 

1^-  SS<M)+  O.S 
SS<M)=  I 
DO  131  M^=  1»  3 

DO  119  1=^  Mf  9»  3 

CALL  EXACT  (SS< I > »PA(M) rCn( I > »B0< I > »N) 

CAI.L  EXACT  (SS(I)  fPB(M>  »CC(I)  »B1(I)  »N) 

CONTINUE 

WRI TF:  <  2 » 1 23 )  <  PA  ( I )  r  I vrj  r  3 )  f  D  r  PB  ( 1 )  » SS  (  1  )  » CC  ( 1  )  » PB  ( 3 )  f  SS (  3  )  r  CC  ( 
H  PB(3>  rSS(3) rCC<3) »D»PP(1  ) fSS(4) fCC 

I  PP(2) »8S(S) tCCCo) rPB(3) »SS(A) »CC(A) 

T  D»PB(1 )fSS(7>rCC<7>tPB(3)rSS(8).rC< 

S  PB(3) fSS(9) »CC(9) 

FORMAT  ('OMFTHOD  DNPl:^ '  » 3< F8 . 4  r  '  N'»7X»'C  ')/ 

B  '  BIN0MIAI..'fFA.lr3<F8.4rF8.1»FB.3)/ 

A  '  ARCSINF  'fFA.lf3(F8.4rF8,l»F8.2)/ 

R  '  NORMAL  '»FA.lr3<F8.4rF8.1rF8.2)) 

WRITE  <1 »13S)  <PA<I)»I-lf3) fD»PB(l) »B0(1 ) i Bl < 1 ) »PB(3) rB0(3> tBI 
P  PB(3)»B0<3)rBl(3)»D»PB(l )»B0(4)»B1 


s  I  Nni.  E » B I  N(5i.  f:/i.  1 : 1.  I  N(ii  p 


F 

T 

S 


0()4?t 

1.25 

F  ORMAT  (  '  OME  THOB  B\P1.  =  '»;?(  F  B .  4  '  P>  1  -A 

M  '  BINOMIAL' »F7.1r9FB*4/'  ARCSINF  ' 

F.  '  NORMAI.  '  f  F  7 . 1 , 9F  8 . 4 ) 

0044 

PB(1.  PB(1)+  riB 

004S 

BPl  -  <riF  0.5)  *PA(2)  -PB<1.  > 

004A 

f; 

00  TO  11* 

0047 

12* 

BF:TA=--  BFTA+  BB 

004B 

IF  (BFTA.I.P  .BH)  00  TO  114 

0050 

BFTA--:  Bl. 

005.1. 

AI.F  A^  AI.FA+  AS 

0052 

IF  (At.FA.I.F  .AB)  00  TO  114 

0054 

URITF  <5»*) 

0055 

URITF  (5»*)'FNri  OF  ITERATIONS!  FNB?  rY/NI 

005* 

Ar.-'Y' 

0057 

RFTURN 

0058 

u 

128 

WRITE  <5f*) 

0059 

write;  (5»*) 'bi-bcrimination  RATIO  error;  n- 

00*0 

A^  'R' 

00*1. 

RETURN 

00*2 

r - 

END 

PB(?)  f»0(S)  »B:i.  (S)  fPBC.?)  »BO<<S)  »B1(A 
nrPBd)  fB0(7)  »BJ.(7)  »PB(;?)  rBO<B)  rBl 
PB(3) »BO(V) »B1 (9) 

P<B  ')/ 


rii 


f  nn 


al4 


n  n  n  n  n  n  n 


siNni.E»siNni.ir;/i.  i :  j  sinoi  k 


0001 


SIJBROUTINF  NORMOI.  (Xt  U»  SOr  7TAB»  TABZ) 


PRORRAM:  NORMAI  probability  DISlRIBimOR  APPROXIMATION 


COM.  NORMAI  -  1  -  F(X)  *  (B1  *  T  +  B:?  ♦  T**:?  +  B3  *  T# 

+  »>!)!(  T#*4  +  Bl=i  #  T# 

MIN.  ACrURACY:  +-  0.000  000  07fi 

RFF:  ABRAMOUITZr  M.»  I.  STFOlJNf  FBS.  'HANriBOOK  OF  MATH. 

FUNCTS.  WITH  FORMOl  AS.  RRAPHS.  ANB  MATH  TAFI.  ES'»  APPI. 
MATH  BFRIFS  tfiS.  WASH.  BC.  NAT'l.  miRKAW  OF  STBS..  19A 


given: 


SAMPl.F  STATIST  in 
=  niSTRIBlITION  MFAN 
-  STANDARD  DEVIATION 


nOMPlJTFS:  7TAB  =  NORMAL.  nUMMlJI.ATIVF  PROBABII.  ITY 
TAB7  ^  <1  -  7TAB> 


BATA  Bl/0.:^19;<sis;?/.  B:?/-0.;iSASA:?7SrV.  B:V1  .7814779S7/.  BlJMY/0 
7.  B4/-l.S:?l?Sf5978/.  B5/1 . ;?;?0?744:?9/ »  GONST/0 . 39894728/ 

7-  <X  -  IJ)/SD 
IF  (7.GE.0.)  00  TO  71 
7.  =  -7 
BUMY  1. 

T  «  l./<l.+7.  ♦0.731A419) 

0  =  1 /SORT <7*  3,141R7A53A>?  SFF  DATA 
F  G*  FXP(  -<Z**2)/7.) 

TAB7.  =  F*  <B1*  T  +B7#  T#<c7  +B3#  7**3  +B4*  T**4  +B5*  T#*S> 

IF  (BUMY.FO.O.)  GO  TO  23 
TAB7  ^  1.-TAP7 

7TAB  «  l.-TABZ 
RETURN 


W  *  V 


a 


F-« 


RIN(5I. 

0001 


P.»P INGLE/LI  :  1=:SINRI.F; 

SUFtROUTTNI:  RFVNOR  (X»  Uf  SP»  7T0»,  TAB7> 


cc— 

nn 

PROGRAM  J  REV 

rr 

rc 

cc 

cc 

MIN 

cc 

cc 

rfe:  hatstini 

cc 

NJr 

C(' 

cc 

given:  1) 

cc 

cn 

cc 

7TAB 

cc 

TAB7 

cc 

cc 

COMPOTFS:  X 

CC-- 

THtl 


:et( 


=  ni SIR I Bin  ION  MEAN 


000 :? 

HAT A  CO/?. 5 1551 7/ » 

Cl /().  802853/ f 

I  Dl/1 .43?788/, 

n2/0.189?A9/» 

ooo:? 

/TAB 

0004 

TOK/-  1-  7TAB 

ooos 

IE  <7tab.gf:.i)  go 

TO  35 

0007 

IE  < /TAB. I  T. 0.5)  GO  TO  31 

ooov 

niJMY=  0. 

0010 

E*  1.-  /TAB 

001 1 

31 

TT=^:  AIOG  (l./E##?) 

001? 

T-=  SORT  <TT) 

0013 

XN“  C0+  Cl*  T+  C?* 

T*»2 

0014 

XB=^  1.+  Bl*  T+  n?* 

T#*?  +B3#  T#« 

0015 

y-  T-  XN/  XB 

001 A 

IE  (BIIMY.EO.O. )  GO 

TO  33 

0018 

X!^^  -X 

0019 

33 

X-  X*  SB4  11 

0070 

GO  TO  37 

00?  1 

35 

Xr-  9V9,999 

00?? 

37 

RETURN 

0023 

C - 

FKB 

al6 


.  tL.'jiil.  Lt  i...  .xiL-Ukrii  'T'“  ":-  *1-  4' 


